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ABSTRACT 

We present high angular resolution Very Large Array multi-epoch continuum observations at 3.6 cm and 7 
mm towards the core of the L723 multipolar outflow revealing a multiple system of four radio sources suspected 
to be YSOs in a region of only ■^4" (1200 AU) in extent. The 3.6 cm observations show that the previously 
detected source VLA 2 contains a close (separation ~ 0"29 or ~90 AU) radio binary, with components (A 
and B) along a position angle of ~ 150°. The northern component (VLA 2 A) of this binary system is also 
detected in the 7 mm observations, with a positive spectral index between 3.6 cm and 7 mm. In addition, the 
source VLA 2A is associated with extended emission along a position angle of ~ 115°, that we interpret as 
outflowing shock-ionized gas that is exciting a system of HH objects with the same position angle. A third, 
weak 3.6 cm source, VLA 2C, that is detected also at 7 mm, is located ^^0^7 northeast of VLA 2A, and is 
possibly associated with the water maser emission in the region. The 7 mm observations reveal the presence of 
an additional source, VLA 2D, located ~3."5 southeast of VLA 2A, and with a 1.35 mm counterpart. All these 
radio continuum sources have a positive spectral index, compatible with them being YSOs. We also propose 
that the high velocity CO emission observed in the region could be the superposition of multiple outflows 
(at least three independent bipolar outflows) excited by the YSOs located at the core, instead of the previous 
interpretations in terms of only one or two outflows. 

Subject headings: ISM: individual (L723) — ISM: jets and outflows — radio continuum: ISM — stars: forma- 
tion 



1. INTRODUCTION 

L723 is an isolated dark cloud at a distance of 300 ±150 
pc (Goldsmith et al. 1984). Located in this cloud is a Class 
source, IRAS 19156H-1906, with a luminosity of -3.4 Lq 
(Dartois et al. 2005). The IRAS source is associated with a 
CO outflow, first mapped by Goldsmith et al. (1984). The 
outflow shows a peculiar quadrupolar morphology, consist- 
ing of two pairs of red-blue lobes with a common center, as is 
clearly seen in the maps of Avery et al. (1990). The larger pair 
of lobes extends along a direction with a position angle (PA.) 
~ 100°, while the smaller pair extends along a direction with 
a PA. ~ 30°. Over the years, the peculiar quadrupolar mor- 
phology of the L723 outflow has been interpreted as due to the 
presence of either one or two driving sources (see Anglada, 
Rodriguez & Torre lies 1996 and references therein). In the 
first case, the quadrupolar morphology is attributed to limb- 
brightening effects in the lobes of a single bipolar outflow, to 
splitting of the lobes due to interaction with ambient clumps, 
or to precession of the outflow axis. In the second case, the 
quadrupolar structure is attributed to two independent bipolar 
outflows, each one driven by a different source. 

Very Large Array (VLA) D-configuration observations (an- 
gular resolution of ^^8") towards the center of the outflow re- 
vealed two 3.6 cm continuum sources, VLA 1 and VLA 2, 
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separated by 15" (Anglada et al. 1991). However, only VLA 
2 is associated with millimeter emission from circumstellar 
dust (Cabrit & Andre 1991), suggesting that it is the counter- 
part of IRAS 19156H-1906. VLA 2 has also been associated 
with a nearby (^O'.'l northeast from VLA 2) H2O maser as 
well as with NH3 emission (Girart et al. 1997). Also, only 
VLA 2 seems to be associated with the CO outflow as shown 
by interferometric observations with the Nobeyama Millime- 
ter Array (Hirano et al. 1998) and with the Berkeley Illinois 
Maryland Association (BIMA) array (Lee et al. 2002). In the 
3.6 cm maps of Anglada et al. (1996), made with data from 
the A and B configurations of the VLA (angular resolution of 
— 0."3), the centimeter emission of VLA 2 is elongated along 
a direction with a PA. ~ 118°, that is close to that of the large 
pair of CO lobes. In addition, H2 emission and Herbig-Haro 
(HH) objects have been found in the region along a direction 
with a PA. ~ 1 15°, similar to that of the elongation of VLA 
2 (Vrba et al. 1986; Hodapp 1994; Palacios & Eiroa 1999; 
Lopez et al. 2006). These results suggest that the outflow ac- 
tivity along the direction ~ 100°-1 18° is driven by the young 
protostellar object VLA 2. However, the identification of the 
exciting source of the smaller pair of CO lobes at P.A.~30° 
remained unclear. As akeady pointed out, one possibility is 
that VLA 2 is the driving source of the two pairs of bipolar 
lobes, e.g., driving each pair in a different epoch as a conse- 
quence of precession of the jet axis. The second possibility is 
the presence of another, still undetected source in the vicinity 
of VLA 2, that could be driving the small pair of CO lobes. 

Other regions harboring multipolar outflows, when studied 
with subarcsecond angular resolution (equivalent to physical 
scales of tens to hundreds of AUs for the relatively nearby re- 
gions that are better studied) and sensitivity at the level of tens 
of fi}y, have shown that the multipolar outflow morphology is 
usually associated with the presence of multiple YSOs at the 
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TABLE 1 
Observations at 3.6 cm 











Bootstrapped 














Effective 


Hux Density of 


Synthesized Beam" 






Observation 


Hour Angle Ranges 


Time on Source 


Phase Calibrator' 


HPBW 


P.A. 


rms Noise* 


Epoch 


Date 


(h) 


(h) 


(Jy) 


(arcsec) 


(deg) 


(/iJy beam"') 


1995.6 


95-Aug-12 


-4A to -2.5 , 2.0 to 4.0 


3.7 


1.146 ± 0.006 


0.32 X 0.27 


-85 


11 


1997.0 


97-Jan-lO 


-2.0 to 2.5 


4.0 


0.926 ± 0.004 


0.27 X 0.25 


15 


11 


1998.2 


98-Mar-27 


1.0 to 3.5, 4.3 to 5.8 


3.7 


0.94 ± 0.01 


0.34 X 0.25 


58 


12 


1998.4 


98-May-26 


1.0 to 3.5, 4.3 to 5.8 


3.7 


0.952 ± 0.007 


0.34 X 0.25 


58 


11 


1999.5 


99-Jul-03 


1.0 to 3.5, 4.3 to 5.8 


3.7 


1.260 ± 0.006 


0.39 X 0.26 


50 


14 


2000.9 


OO-Dec-16 


-5.0 to 4.5 


8.0 


1.65 ± 0.02 


0.28 X 0.27 


-70 


6 



^The phase calibrator used in all the observations was 1925+211. 
''For naturally weighted maps. 



TABLE 2 
Observations at 7 mm 







Bootstrapped 
Flux Density of 


Synthesized Beam*" 






Observation 


VLA 


Phase Calibrator" 


HPBW 


RA. 


rms Noise'' 


Date" 


Configuration 


(Jy) 


(arcsec) 


(deg) 


(mjy beam"') 


02-Jan-05 


D 


2.38 ± 0.04 


2.37x1.75 


-22 


0.12 


04-Feb-Ol 


CnB 


1.85 ±0.05 


0.45x0.39 


-76 


0.12 


04-Mar-26 


C 


2.22 ± 0.05 


0.58x0.51 


-8 


0.12 


05-Jul-29 


C 


1.68 ± 0.02 


0.57x0.51 


+34 


0.12 



"The phase calibrator used in all the observations was 1925+21 1 . 
''For naturally weighted maps. 



center of the outflows: e.g., IRAS 16293-2422 (CO outflow: 
Mizuno et al. 1990; radio continuum: Estalella et al. 1991); 
IRAS 20050+2720 (CO outflow: Bachiller et al. 1995; radio 
continuum: Anglada, Rodriguez & Torrelles 1998); HHlll 
(CO outflow: Cernicharo & Reipurth 1996; radio continuum: 
Reipurth et al. 1999); L1634 (H2 jet: Hodapp & Ladd 1995; 
radio continuum: Beltran et al. 2002); HH288 (CO outflow: 
Gueth et al. 2001; radio continuum: Franco-Hernandez & 
Rodriguez 2003); OMC-IS (optical jets: BaUy et al. 2000; 
radio continuum: Zapata et al. 2004). 

In this paper we present sensitive, high angular resolution 
multi-epoch observations at 3.6 cm of the radio sources de- 
tected in the central region of the L723 outflow. Our goals 
were to search for changes in the structure of VLA 2 that 
could help us to clarify the origin of the quadrupolar mor- 
phology of the outflow and to search for additional indica- 
tions of multiplicity at very small scales. We also analyze 
VLA archive 7 mm continuum data of the region. 

2. OBSERVATIONS 

The 3.6 cm continuum observations were carried out in six 
epochs, ranging from 1995 August 12 to 2000 December 16. 
The data from the 1995.6 observations were already reported 
by Anglada et al. (1996). The VLA of the National Radio 
Astronomy Observatory (NRAO)^ was used in its A config- 
uration, providing an angular resolution of ^0"3. Absolute 
flux calibration was achieved by observing 3C 286, for which 
a flux density of 5.20 Jy at 3.6 cm was adopted. Data edit- 
ing and calibration were carried out using the Astronomical 
Image Processing System (AIPS) package of NRAO, follow- 
ing the standard VLA procedures. The observation dates, the 

^ The NRAO is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 



bootstrapped flux densities of 1925+211, the phase calibra- 
tor used in all the observations, as well as the parameters 
of the synthesized beams, and the rms noise of the naturally 
weighted maps are given in Table 1 . Due to the weakness of 
the sources detected the signal-to-noise ratio per baseline was 
insufficient to self-calibrate the data. 

For each epoch, we made cleaned maps with different val- 
ues of the parameter ROBUST of IMAGR (Briggs 1995). 
The ROBUST parameter ranges in value from +5 to -5, 
with the maps made with ROBUST=+5 (equivalent to natural 
weighting of the visibiUty data) having the highest sensitiv- 
ity but the lowest angular resolution and the maps made with 
ROBUST=-5 (equivalent to uniform weighting) having the 
highest angular resolution but the lowest sensitivity. For each 
value of the ROBUST parameter the map of each epoch was 
restored using a circular beam whose HPBW is the geometric 
mean of the major and minor axes of the synthesized beam 
obtained by concatenating the uv data of all the six epochs 
(HPBW = 0."29 for naturally weighted maps and HPBW = 
Of 19 for uniformly weighted maps). This average restoring 
beam was adopted to facilitate comparison between the differ- 
ent epochs and to allow averaging of all the 3.6 cm images into 
a single image. We consider this procedure appropriate since 
the major and minor axes of the synthesized beams of the in- 
dividual observing sessions (see Table 1) are within ~10% of 
the adopted average HPBW. 

To improve the sensitivity and uv coverage, we have con- 
catenated the uv data from all the epochs to obtain single 
maps, which we will refer to as "All-Epochs" maps. Since 
some of the structures persist from one epoch to another, 
while others change, the procedure of concatenating data 
could result in spurious effects in the final images. There- 
fore, we checked that the maps obtained by concatenating all 
the uv data do not show significant differences from the maps 
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TABLE 3 

YSO Candidates Associated with VLA 2 



Component 




Position (12000)" 




Flux Density (mjy)" 


Spectral 
Index 


RA 


DEC 


Error 


3.6 cm 


7 mm 


A 


19 17 53.673 


+19 12 19.59 


o;'oi 


0.11 ±0.04'= 


0.7 ± 0.2 


+1.1 ±0.3 


B 


19 17 53.685 


+19 12 19.34 


O'.'Ol 


0.08 ± 0.02 


<0.37'' 


<+0.9 


C 


19 17 53.739 


+19 12 19.77 


O'.'Ol 


0.05 ± 0.02= 


0.7 ± 0.2f 


+1.5 ±0.4 


D 


19 17 53.894 


+19 12 17.96 


0"09 


<0.02'' 


0.6 ± 0.2 


>+2.0 



"Positions for components A and B are derived from Gaussian fits to the sources in the 3.6 cm "AU-Epochs" imiformly weighted map (Fig. 3a). Positions for components C and D 
are derived from the 7 mm map (Fig. 4). Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. 
''Derived from Gaussian fits to "All-Epochs" naturally weighted maps. 

'= Uncertainty is estimated from the observed pealc intensity variabihty in our multi-epoch observations (Table 4) 

upper limit. 

■^Sum of the flux densities of VLA 2Ca and VLA 2Ch. 

^Angularly resolved source. Deconvolved size is ~ 1^'2 ± at a PA 2:; 90° ±7°. 



obtained by averaging all the single epoch maps. 

The 7 mm continuum data were taken from the VLA 
archive. The data consist of four sets of observations in the 
D, C and CnB configurations (see Table 2). Absolute flux 
calibration was achieved by observing 3C 286, and a model 
image at 7 mm for this calibrator was used to cahbrate the data 
in amplitude. In order to improve the sensitivity of the final 
image, but keeping a balanced contribution from the different 
configurations, we have concatenated all the available uv data 
giving the same weight to the D and C+CnB data. Since these 
data come from different configurations, it is not possible to 
average the individual maps. From these uv data we obtained 
a 7 rmn "AU-Epochs" naturally weighted map. 

3. RESULTS AND DISCUSSION 

3.1. VLA 2A and2B: A Close Binary System ofYSOs 

In Figure 1 we show naturally weighted 3.6 cm maps of 
VLA 2 for each individual epoch. We have detected radio 
continuum emission from VLA 2 in all the epochs, with total 
flux densities varying in the range 0.2-0.5 mJy. In these maps, 
we see changes in the morphology of VLA 2, but some of 
the structures persist from one epoch to another. In particular, 
in three of the epochs VLA 2 shows a clear double morphol- 
ogy with two components separated by ~0"3 along a P. A. 
of ~150°. This double morphology appears in the 1998.2 
and 1998.4 epochs, where the uv coverage was exactly the 
same, and in the 2000.9 epoch, which was a full-track ob- 
servation and therefore the uv coverage obtained was much 
better than in the other epochs (see Table 1 for the HA ranges 
and time on source for each epoch of observation). In the 
1999.5 epoch, the uv coverage was also very similar to that 
of the 1998 epochs, but the double morphology is not so evi- 
dent. We think this could be due to the higher rms noise of the 
1999.5 map and to the increase of the flux density of the north- 
ern component (see below) that partially masks the emission 
of the southern component. Given these considerations, it is 
likely that VLA 2 is actually a persistent double radio source, 
with the differences in the observed morphology being mainly 
due to the different uv coverage from one epoch to another. 

A way to mitigate the effects introduced by the different uv 
coverages would be to make uniformly weighted maps of the 
individual epochs. In uniformly weighted maps all the grid 
cells have equal weight, which compensates for the high con- 
centration of data points in the inner region of the uv plane. 
This provides a more uniform coverage of the gridded uv 
plane, resulting in maps having higher angular resolution than 
naturally weighted maps. Unfortunately, uniform weighting 



TABLE 4 

Monitoring of the VLA 2 binary at 3.6 cm 



Epoch 


Peak Intensity (mJy beam ' )" 


Total Flux" 
Density (miy) 


VLA 2A 


VLA2B 


1995.6 


0.08 ± 0.01'^ 


0.06 ± 0.01'^ 


0.30 ± 0.03 


1997.0 


0.07 ± 0.01'^ 


0.06 ± 0.01'^ 


0.18 ± 0.03 


1998.2 


0.06 ± 0.01 


0.06 ± 0.01 


0.20 ± 0.02 


1998.4 


0.06 ± 0.01 


0.07 ± 0.01 


0.16 ±0.02 


1999.5 


0.15 ± O.OP 


0.05 ± O.OP 


0.21 ± 0.03 


2000.9 


0.13 ±0.01 


0.07 ± 0.01 


0.24 ± 0.02 



"Derived from Gaussian fits to the 3.6 cm naturally weighted maps that have been 
restored with an average circular beam of HPBW=()''29 (Fig. 1). 

'^Totai flux density of the binary system excluding the extended emission (see Fig. 1). 

""Obtained by fixing the Gaussian centers to the positions of VLA 2A and VLA 2B 
given in Table 3. 



yields lower sensitivity, and due to the weakness of the source, 
it is difficult to detect it in the uniformly weighted map of a 
single epoch. However, it is still possible to test our hypoth- 
esis of a double radio source with the available data. For this 
purpose we have made uniformly weighted maps from two 
different sets of concatenated uv data: (i) the three epochs 
in which VLA 2 does not show clearly a double morphology, 
and (ii) the three epochs in which VLA 2 shows a clear double 
morphology. This procedure allows us to compare two maps 
with high enough angular resolution and sensitivity, as well 
as a good uv coverage. As can be seen in Figure 2, VLA 2 
shows always a double morphology (with a slight change of 
~0"1 in the position of the sources), even in the map made 
with the data of the epochs in which the double morphology 
was not evident. This strongly supports that two radio con- 
tinuum sources (that we call VLA 2A and VLA 2B) persist 
through all the epochs, covering a time span > 5 yr. 

Since all the data seem to be consistent with the presence 
of a double source, we obtained an "All-Epochs" uniformly 
weighted map (Fig. 3a) by concatenating all the uv data. The 
positions and flux densities of components VLA 2A and 2B, 
derived from Gaussian ellipsoid fits to this map, are given in 
Table 3. The angular separation of the two components is 
0"29, which corresponds to a projected separation of ^90 
AU (assuming a distance of 300 pc). If these sources were 
tracing some kind of high-velocity outflowing features, such 
as knots in a jet, a significant variation in their positions in 
the single epoch maps would be expected (e.g., a velocity of 
'--^100 km s"' in the plane of the sky would result in an angular 
displacement of ~ 0"35 over the period of ~5 yr covered by 
our observations). The lack of a significant displacement from 
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J I I L 



1998.4 
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(d) 
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n I I I r 



1999^5^ 




(e) 




(f) 



1917 53.78 53.76 53.74 
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Fig. 1. — VLA 3.6 cm continuum maps of the VLA 2 region at different epochs, obtained with natural weighting. In order to better identify morphological 
changes from one epoch to other, the maps have been restored using the same circular beam of Of' 29 (HPBW) and contoured in terms of the same reference value 
(11 /^Jy beam~^). For panels (a) to (e), contour levels are —4, —3, 3, 4, 5, 6, 7, 8, 9, 10, and 12 times 1 1 /jJy beam~^ the typical rms of these maps (see Table 
1). For panel (f), additional contours at -2, —1.5, L5,and2times 11 /^Jybeam"^ have been included so that the lowest contour level corresponds approximately 
to 3 times the rms of this map (see Table 1). The crosses mark the positions of VLA 2A and VLA 2B derived from Gaussian fits to the uniformly weighted 
'All-epochs" map (see text). The filled square marks the position of file H2O maser detected by Girart et al. (1997). The radio continuum knots detected in the 
region (w2, wl, el, e2, si, and nl) are also labeled (see §3.1). 
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L723 VLA 2 
3.6 cm 



I I r 

1995.6+1997.0+1999.5 



(a) 



18.4 I 



O 



I I r 

1998.2+1998.4+2000.9 




(b) 



Fig. 2. — (a) Unifomily weighted map of the source VLA 2 at 3.6 cm 
made by concatenating the uv data of the 1995.6, 1997.0, and 1999.5 epochs 
(where VLA 2 does not show a clear double morphology; see Fig. 1). The 
synthesized beam of the map is 0^'20xOy 18, with a P.A. of 79°. Contour 
levels are -4, -3, 3, 4, and 5 times the rms of the map, 17 ^l}y beam"^(b) 
Same as (a), but for the 1998.2, 1998.4, and 2000.9 epochs (where VLA 2 
shows a clear double morphology; see Fig. 1). The synthesized beam of the 
map is 0"20xO"18, with a P.A. of 86°. Contour levels are -4, -3, 3, 4, 
and 5 times the rms of the map, 12 /ijy beam"' . Note that in both maps the 
source VLA 2 shows a clear double morphology, suggesting that VLA 2 is a 
binary system of embedded YSOs (see text for further details). The crosses 
mark the positions of VLA 2A and VLA 2B derived from Gaussian fits to the 
unifomily weighted "All-epochs" map (see text, and Fig. 3a). 

one epoch to another excludes that these sources trace high- 
velocity material, unless the axis of the motion was very close 
the Une of sight (which does not seem to be the case, given 
that the blue and red lobes of the molecular outflow appear 
clearly separated on the sky). This, together with the addi- 
tional properties observed for these sources (see below) led 
us to interpret VLA 2A and VLA 2B as a close binary system 
of embedded YSOs. 

The individual epoch maps shown in Figure 1 suggest vari- 
ability in the 3.6 cm flux density of the compact components 
of the double source. In Table 4 we give the peak intensities 
of VLA 2A and VLA 2B at each epoch of observation. Be- 
cause the angular resolution of our observations is similar to 
the angular separation of the two sources, these values of the 
peak intensity must be taken with care. However, given that 
the maps of the different epochs have been restored with the 



same beam, the variations of the peak intensities are expected 
to be representative of the variability of the individual compo- 
nents. The peak intensity of VLA 2B appears to be constant at 
aU epochs, consistent with a value of 0.06 ± 0.01 mJy beam~^ 
(see Table 4). On the other hand, the peak intensity of VLA 
2A is consistent with a value of 0.07 ± 0.01 mJy beam"' (see 
Table 4) in the four first epochs, while in the last two epochs 
the peak intensity of the source increased approximately by 
a factor of two. Similar flux density variability has been ob- 
served in several embedded YSOs (e.g., Cep A, Porras et al. 
2002; B335, Avila et al. 2001). We interpret the variabil- 
ity in the flux density of VLA 2A as due to material recently 
ejected that still remains very close to the YSO and it is not 
distinguishable from the source. 

In addition, we have detected 7 mm emission from VLA 2A 
(see Fig. 4). The flux densities at 3.6 cm and 7 mm of VLA 
2A (see Table 3) imply a spectral index a ~ 1.1 ±0.3 (S,^ oc 
p"), consistent with partially optically thick free-free emis- 
sion that could be tracing a radio jet. The nominal value of 
the spectral index obtained is slightly higher than the typical 
values found in ionized jets driven by YSOs, 0.3 < a < 0.9 
(Anglada 1996 and references therein). This higher value of 
the spectral index of VLA 2A could be indicating that we are 
partially missing the extended emission of the possible radio 
jet, either because of insufficient uv coverage or because part 
of this extended emission is confused with the emission of the 
nearby source VLA 2B. Since the outer part of a thermal ra- 
dio jet is where the emission is optically thinner, we would 
be observing preferentially the emission coming from the re- 
gion closer to the driving source, where the optical depth is 
higher. This effect would increase the value of the spectral 
index with respect to that of the overall emission of a radio 
jet. Alternatively, a high value of the spectral index could be 
indicative of a significant contribution from dust emission at 
7 mm. Both possibilities (free-free emission from a jet with 
a possible contribution of dust emission from a circumstel- 
lar disk or an envelope) are consistent with VLA 2A being a 
YSO. 

We should note that the spectral indices derived in this pa- 
per have been obtained from non-simultaneous observations 
and can be affected by time-variability of the sources. How- 
ever, as discussed above, the flux variations observed in our 
multi-epoch observations suggest modest variability that it is 
not expected to affect our conclusions on the nature of the 
emission derived from the spectral index. Also, the 3.6 cm 
and 7 mm observations have different uv coverages, that could 
introduce additional errors in the flux density measurements 
of extended emission. Finally, the presence of several sources 
in a small area makes difficult to isolate the contribution of 
each source. Therefore, the spectral indices reported in this 
paper should be considered as rough estimates. 

In addition to the compact emission of the two sources VLA 
2 A and VLA 2B, in the individual 3.6 cm maps of Figure 1 
we also note the presence of several weak "knots" that could 
be ttacing extended emission. In order to gain sensitivity to 
extended structures, we have concatenated all the 3.6 cm uv 
data and made a naturally weighted "All-Epochs" map, that 
is shown in Figure 3b. This map shows extended emission 
with a P.A. of ~115°, where we identify the knots wl, w2, 
el, and e2, that can also be seen in several of the individual 
epoch maps (Fig. 1). The extended emission seems to em- 
anate from VLA 2A, the northern component of the double 
source, since a line running through the major axis of this ex- 
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Fig. 3. — (a) Uniformly weighted "All-Epoclis" map of the source VLA 2 at 3.6 cm, obtained by concatenating the uv data from all the epochs. The synthesized 
beam of the map is ()!'2()x()" 18, with a P.A. of 86°. Contour levels are -4, -3, 3, 4, 5, and 6 times the rms of the map, 1 1 /iJy beam"' . (b) Same as (a), but with 
natural weighting. The synthesized beam of the map is 0;'31 x()"27, with a P.A. of 63°. Contour levels are -4, -3, 3,4, 5, 6, 7, 8, 9, 10, 12, 14, 16, 18, and 20 
times the rms of the map, 4 fi]y beam"' . Crosses mark the positions of VLA 2A and VLA 2B derived from Gaussian fits to the uniformly weighted map. The 
filled square marks the position of the H2O maser detected by Grrart et al. (1997). The axis of the extended emission (traced by knots w2, wl, el, and e2), at 
P.A. of ~115°, is indicated. 



tended emission is centered close to the position of VLA 2A 
(see Fig. 3b), which is the strongest source along this direc- 
tion. The P.A. of this extended emission (^115°) coincides 
with that of the system of HH objects observed in the region 
(Lopez et al. 2006), suggesting that VLA 2A is a YSO that is 
driving these HH objects. The P.A. of the extended emission 
is also similar to that of the symmetry axis of the larger pair 
of CO outflow lobes (-100°), thus suggesting that VLA 2A 
is also associated with this lobe pair. 

We interpret the extended radio emission as tracing an ion- 
ized jet that originates at the position of VLA 2A, and where 



the knots wl, w2, el, and e2 would be tracing electron den- 
sity enhancements produced by shock-ionization. We do not 
detect clear proper motions associated with these knots from 
one epoch to another, setting an upper limit to the proper mo- 
tion velocity of the knots of ~27 km s"^ (however, due to 
the coarse temporal sampling and the low signal-to-noise ra- 
tio of the knots in the individual epoch maps, they are diffi- 
cult to identify from one epoch to another, and we cannot dis- 
card higher proper motion velocities). This suggests that these 
knots have been excited by the interaction of the jet of VLA 
2A with dense, stationary clumps in the surrounding gas. 
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We also note the detection of the faint knots, nl (in the 
2000.9 image; Fig. If) and si (in the 1997.0 image; Fig. 
lb). We interpret these two knots as shock-ionized gas as- 
sociated with episodic ejections of material. Their location, 
nearly aligned with source VLA 2B, at a P.A. ~ 20°, is con- 
sistent with being excited by a jet coming from this source. 
The P.A. of this jet is similar to that of the smaller pair of CO 
outflow lobes (~30°) and therefore, we tentatively propose 
VLA 2B as its exciting source. 

In summary, we detected two compact 3.6 cm sources, 
VLA 2A and VLA 2B, separated by ~0."3 along a P.A. of 
~150°. We also detected extended 3.6 cm emission, centered 
on VLA 2A, and along a P.A. of ~115°, that coincides with 
the P.A of the system of HH objects in the region and is simi- 
lar to that of the larger pair of CO outflow lobes. The spectral 
index of VLA 2 A in the 3.6 cm to 7 mm wavelength range, is 
consistent with this source tracing the origin of an ionized ra- 
dio jet with a possible contribution from dust. Therefore, we 
conclude that VLA 2A is an embedded YSO associated with 
a radio jet; it is also the likely exciting source of the system 
of HH objects and of the larger pair of lobes of the CO out- 
flow observed in the region. We identify several knots (wl, 
w2, el, and e2) distributed along the extended 3.6 cm emis- 
sion, at both sides of VLA 2A. Since these knots are much 
weaker than VLA 2A and are aligned along the axis of the 
extended emission we interpret them as tracing electron den- 
sity enhancements in the jet from VLA 2 A rather than tracing 
the location of other YSOs. The lack of clear proper motions 
suggests that these knots have been excited by the interaction 
of the jet with dense, stationary clumps in the gas surround- 
ing VLA 2A. The compact source VLA 2B has a flux density 
comparable to that of VLA 2A and both sources are ahgned 



along a P.A. of ~150°, a direction where there are no indi- 
cations of outflow activity. Therefore, we interpret VLA 2B 
as a second embedded YSO, constituting a close (projected 
separation of 90 AU) binary system. Neither VLA 2A nor 
VLA 2B show noticeable proper motions; this result strongly 
suggests that these sources do not trace high-velocity mate- 
rial, and is consistent with these sources being YSOs. Un- 
fortunately, with the current data it is not possible to obtain 
a precise enough measure of proper motions that would al- 
lows us to derive possible orbital motions, as has been done 
in other sources (e.g., Loinard et al. 2002); this kind of data 
could be helpful to confirm that this is a binary system. Also, 
there are no high angular resolution data outside the centime- 
ter wavelength range that allow us to characterize the individ- 
ual spectral energy distribution of each component in order to 
establish the evolutionary status of the sources and to confirm 
their nature as YSOs. Therefore, these sources are good tar- 
gets for future observations with the EVLA and ALMA. We 
also detect in the single epoch images two faint knots (nl and 
si) nearly aligned with the source VLA 2B along a P.A. of 
^20°, that is similar to the P.A. of the smaller pair of CO out- 
flow lobes. We tentatively interpret these two knots as shock- 
ionized gas associated with episodic ejections from VLA 2B, 
that could be the driving source of these knots as well as of 
this pair of CO outflow lobes. 

3.2. Other Radio Sources Detected in the Region 

In addition to VLA 2A, our 7 mm image (Fig. 4) shows 
emission from other two sources. One of them, that we will 
refer to as VLA 2C, coincides with the position of the water 
maser emission detected by Girart et al. (1997) and Furuya 
et al. (2003). The coincidence of this source with the water 
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TABLE 5 
Monitoring of VLA 1 at 3.6 cm^ 



Deconvolved Size 







Position (J2000)'' 




Hux Density 


Major Axis 


Minor Axis 


P.A. 


Epoch 


RA 


DEC 


Error 


(mjy) 


(arcsec) 


(arcsec) 


(deg) 


1995.6 


19 17 52.922 


+19 12 08.85 


o;'oi 


0.36 ± 0.02 


0.11 


0.10 


179 


1997.0 


19 17 52.923 


+19 12 08.85 


O'.'Ol 


0.38 ± 0.02 


0.13 


0.00 


37 


1998.2 


19 17 52.922 


+19 12 08.86 


O'.'Ol 


0.34 ± 0.02 


0.12 


0.04 


179 


1998.4 


19 17 52.922 


+19 12 08.86 


O'.'Ol 


0.31 ± 0.02 


0.11 


0.05 


179 


1999.5 


19 17 52.923 


+19 12 08.86 


O'.'Ol 


0.34 ± 0.03 


0.13 


0.07 


179 


2000.9 


19 17 52.923 


+19 12 08.85 


O'.'Ol 


0.34 ± 0.01 


0.13 


0.00 


30 



^Parameters derived from Gaussian fits to the 3.6 cm naturally weighted maps that have been restored with an average circular beam of HPBW=0f'29. 
'^Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcininutes, and arcseconds. 



maser suggests that VLA 2C is a YSO. At 7 mm, the source 
VLA 2C appears elongated, with a deconvolved size of 1 "2 
± O'.'l with a P.A. ~ 90° ± 7°. At 3.6 cm, we detect weak ex- 
tended emission that coincides with the source detected at 7 
mm. The 3.6 cm emission shows two peaks, that we call VLA 
2Ca and VLA 2Cb. With the current data it is unclear if these 
two peaks actually correspond to two objects or to a single ex- 
tended source. Sensitive high angular resolution multiwave- 
length observations would be required to further define the 
nature of this emission. From the total flux density at 3.6 cm 
(VLA 2Ca + VLA 2Cb) and 7 mm (see Table 3), we estimate 
a spectral index of 1.5 ±0.4 for VLA 2C, a value suggestive 
of thermal emission from a YSO (see discussion in §3.1). 

The third 7 mm source, VLA 2D, is located ^ 3 "5 from 
VLA 2 A and does not have a detected 3.6 cm counterpart. 
However, this source is located within ^0"6 of the nominal 
position of a 1.35 mm source detected by Girart et al. (2007) 
with the SubmilUmeter Array (SMA) (see Fig. 4). Given 
that the beam size of the SMA observations is 3"3xl."5 
(P.A. =79°), we suggest that VLA 2D is the 7 mm counterpart 
of the 1.35 mm source. From the total flux density at 7 mm 
of VLA 2D, and adopting a 4a upper limit to the 3.6 cm flux 
density, we estimate a lower Umit to the spectral index of this 
source of ^2.0 in this wavelength range, suggesting that the 
7 mm emission is dominated by thermal dust emission. Ad- 
ditionally, the 1.35 mm source seems to be associated with an 
elongated SiO structure, probably tracing material shocked by 
a high velocity flow (Girart et al. 2007), suggesting that this 
source is associated with a YSO. 

In addition to the sources VLA 2A, VLA 2B, VLA 2C, and 
VLA 2D that appear located at the center of the L723 outflow, 
our 3.6 cm images detect the source VLA 1, first reported by 
Anglada et al. (1991), and located ~15" SW of the position 
of VLA 2A. The source VLA 1 appears always compact (de- 
convolved size < 0"13; see Table 5) with a roughly constant 
flux density consistent at all epochs with a value of 0.36±0.02 
mJy (see Table 5). We have not detected this source at 7 mm 
up to a 4cr level of 0.32 mJy. With our 3.6 cm flux density 
and the 7 mm upper Umit, we estimate for VLA 1 a spec- 
tral index a < -0. 1 , a value consistent with the spectral index 
q;=-0.3±0.2 measured by Anglada et al. (1996) between 3.6 
and 6 cm. As discussed by Anglada et al. (1996) and Girart 
et al. (1997), this source is most probably a radio-emitting, 
optically-obscured T Tauri star located very near to the denser 
part of the molecular cloud. Our accurate astrometry indicates 
that this compact source does not show noticeable proper mo- 
tions, with a position displacement < 0."03 (3f7) over the 5.3 
yr period covered by our observations (see Table 5). Since 
a total proper motion displacement of ~ 0!'03 is expected for 



an object embedded in the L723 cloud (assuming a distance of 
300 pc and the Brand & Blitz 1993 model of galactic rotation), 
this result is compatible with VLA 1 belonging to the L723 
cloud. However, the possibility that VLA 1 is a line-of-sight 
source not directly associated with the molecular cloud is not 
discarded. The remaining radio sources (VLA 2A, VLA 2B, 
VLA 2C, and VLA 2D) discussed in this paper are embedded 
in regions of extended radio emission and are associated with 
ejecta, which does not allow a precise determination of their 
possible proper motions. 

3.3. The Morphology of the Outflow in L723: Three Bipolar 

Outflows? 

In the maps of Avery et al. (1990), the L723 CO outflow 
shows a four-lobe morphology with a large pair of bipolar 
lobes along a P.A. of ~ 100° and a smaller pair of bipolar 
lobes along a P.A. of ^ 30°. As discussed in the previous sec- 
tions, our radio continuum data suggest that a multiple system 
of at least four YSOs is present, within a region of ~ 1200 AU, 
at the intersection of the axes of symmetry of the two pairs of 
CO lobes. Furthermore, in §3.1 we discussed the results that 
led us to propose that VLA 2A is related to the excitation of 
the larger pair of CO lobes, while VLA 2B could be driving 
the smaller pair of CO lobes. 

Nevertheless, the interferometric maps of Lee et al. (2002) 
(obtained with higher angular resolution than the maps of Av- 
ery et al. 1990) reveal additional details of the CO outflow 
structure (see Fig. 5a). As can be seen in this figure, the 
smaller pair of CO lobes shows a compact shape with a well 
defined axis at P.A. ~ 30°. On the other hand, the shape of the 
larger pair of CO lobes appears clumpy, with only a few weak 
CO emission knots near the symmetry axis of this lobe pair 
(at P.A. ~ 100°), and most of the emission concentrated at 
both sides of this central axis, roughly at P.A. ~ 90° and PA. 

115°. The distribution of the CO emission of this larger 
lobe pair was interpreted by Lee et al. (2002) as tracing the 
Umb-brightened cavity walls of a bipolar outflow with a P.A. 
~ 100°. However, we note that the distribution of the CO 
emission at both sides of the central axis (P.A. ~ 100°) of this 
lobe pair is asymmetric. The emission at PA. ~ 115° starts 
at the center of the outflow (the position of VLA 2) and is 
stronger than the emission at P.A. ~ 90°, which shows a gap, 
with the bulk of the emission apparently starting ~ 1 ' away 
from the outflow center Furthermore, the P.A. of 1 15° coin- 
cides with that of the extended radio emission associated with 
VLA 2A and with that of the system of HH objects. 

These considerations led us to propose a new interpreta- 
tion of the high-velocity CO emission in L723 in terms of 
three bipolar outflows. One bipolar outflow at P.A. ~ 30°, 
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Fig. 5. — {a) Map of the high velocity CO (/=!— >0) emission in L723 (Lee et al. 2002). Redshifted emission (dashed contours) is integrated from 11 to 31 
km s~' and blueshifted emission (solid contours) is integrated from —10 to 10 km s"'. Dashed lines mark the position angles of the axes of the bipolar outflows 
proposed by Lee et al. (2002). Sohd lines mark the position angles of the axes of the three bipolar outflows proposed in this paper. Herbig-Haro objects (Lopez 
et al. 2006) are marked by boxes. The star symbol marks the position of VLA 2 and the triangle mai'ks the position of VLA 1 (Anglada et al. 1991, 1996). (b) 
A schematic depiction of the proposed outflow axes and the multiple system of YSOs at the center of the L723 outflow. Stai' symbols mark the positions of the 
suspected YSOs: VLA 2A, VLA 2B, VLA 2C, and VLA 2D. The thi'ee axes have the directions of the outflows proposed in this paper. The crosses mark the 
positions of the knots detected at 3.6 cm. 
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that could be driven by VLA 2B. A second bipolar outflow at 
P.A. ~ 1 15°, which is associated with a system of HH objects 
and shock-ionized extended radio continuum emission with 
the same position angle, that is driven by VLA 2A. Finally, a 
third bipolar outflow at P.A. ~ 90°, whose exciting source is 
uncertain. The lobes of the proposed outflows at P.A. ~ 30° 
and P.A. ^115° appear well "connected" with the outflow 
center, suggesting recent outflow activity in these directions. 
On the contrary, the lobes of the outflow at P.A. ~ 90° are 
weaker and weU detached from the central region, suggest- 
ing that they resulted from activity taking place in the past. 
Therefore, we speculate that this outflow at P.A. ~ 90° could 
be a "fossil" outflow in the sense that its exciting source has 
not been active in the recent past. The motion oljserved will 
then be the result of momentum conservation of the gas accel- 
erated when the exciting source was active. A similar fossil 
molecular outflow has been proposed to exist in the molecular 
cloud B5 (Yu, BiUawala, & Bafly 1999). 

4. CONCLUSIONS 

Our high angular resolution 3.6 cm and 7 mm continuum 

observations towards the previously known source VLA 2, at 
the center of the L723 multipolar CO outflow, have revealed a 
group of at least four YSOs within a region of ~ 4" (1200 AU) 
in size. The size of this region is comparable with the size of 
the region where Girart et al. (1997) found a local heating of 
the high density gas. Two of our continuum sources, VLA 2A 
and VLA 2B, form a close radio binary whose components are 



separated by --0."29 90 AU). Another source, VLA 2C, is 
associated with a water maser, and a fourth source, VLA 2D, 
that we only detect at 7 mm, appears to be the counterpart of 
a recently discovered 1.35 mm source. 

We propose that the multipolar CO outflow in L723 could 
result from the superposition of at least three independent 
bipolar outflows with position angles of ^115°, ^90°, and 
~30°, driven by three different YSOs. Our 3.6 cm observa- 
tions suggest that VLA 2A is associated with an ionized jet at 
a position angle of ^ 1 1 5 ° . We propose that this is the exciting 
source of the system of HH objects previously detected in the 
region, as well as of a bipolar CO outflow with the same po- 
sition angle. The presence of two 3.6 cm knots aligned with 
VLA 2B at a P.A. of --20° suggests that VLA 2B could be the 
exciting source of the bipolar CO outflow at a P.A of ~30°. 
Finally, the third bipolar CO outflow, with a position angle of 
~ 90°, seems to be a "fossil" outflow whose exciting source 
has not been very active in the recent past. 
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